Do attractive bosons condense? 
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Motivated by experiments on bose atoms in traps which have attractive interactions (e.g. '^Li), 
we consider two models which may be solved exactly. We construct the ground states subject to the 
constraint that the system is rotating with angular momentum proportional to the number of atoms. 
In a conventional system this would lead to quantised vortices; here, for attractive interactions, we 
find that the angular momentum is absorbed by the centre of mass motion. Moreover, the state is 
uncondensed and is an example of a 'fragmented' condensate discussed by Nozieres and Saint James. 
The same models with repulsive interactions are fully condensed in the thermodynamic limit. 



One of the most novel aspects of the creation of Bose 
condensates with neutral atoms in traps is the possibil- 
ity of observing a bose gas with attractive interactions 
(negative scattering lengths). The case of ^Li has been 
studied both experimentally and theoretically. Con- 
densation has been predicted to be stable for a sufficiently 
small number of particles or sufficiently weak interactions 
The instability to collapse when these conditions 
are not obeyed has also been discussed by several authors 

In this Letter we show, using two exactly soluble 
models, that there may be other possibilities for non- 
condensed states with attractive interactions. The states 
are the 'fragmented' condensates discussed by Nozieres 
and Saint James [ |To| in the context of excitonic bose con- 
densates. The possibility of such states emerges from the 
realisation that it is the exchange interaction which 
causes bosons with repulsive interactions to condense into 
a single one-particle state, if there are several one-particle 
ground states. Conversely for attractive interactions, the 
exchange term is negative and may prefer 'fragmented' 
pof condensation into more than one state if there is a 
degeneracy (or perhaps if the interactions are sufficiently 
strong). Kagan et al. [|| argue that trapped gases with 
sufficiently large negative scattering lengths are unstable 
to the formation of clusters using a somewhat different 
argument, but with the same physical origin. 

The two models we examine are: particles in a har- 
monic trap with L quanta of angular momenta with at- 
tractive interactions treated as a degenerate perturba- 
tion [|2j; rotating particles in a harmonic trap interact- 
ing with harmonic interactions (Both of these 
cases have been of interest for fermions |12jl3|, where 



rotation is replaced by a magnetic field and the phenom- 
ena are related to the fractional quantum Hall effect.) 
Rotation is considered in both cases, partly because the 
non-rotating ground state, in the thermodynamic limit, 
is trivial in both cases (for different reasons) and partly 



because the response to rotation is characteristic of su- 
perfluidity in the system . 

Consider the two-dimensional Hamiltonian 
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in the limit where the dimensionless coupling is weak, 
1 77 1 <C 1, so that the contact interaction can be treated 
perturbatively. We will now determine the ground state 
subject to the constraint that the system contains L 
quanta of angular momentum. We should note that the 
centre of mass variables will separate in this Hamiltonian 
because the trap is harmonic. This will be used below. 

The single particle spectrum is usefully expressed |l9| ] 
in terms of the angular momentum quantum number, m, 
and the radial quantum number, n,-: 



E = |m| + 2nr + 1 



(2) 



in dimensionless units. For the non-interacting case, 
7y = 0, with N particles it is clear that to minimise the 
energy n-^ = Q for all particles. The angular momentum 
may be expressed as L = mNm where Nm is the 
number of particles in the state m. There is a degen- 
eracy, in general, associated with the choice of the set 
{Nra\- If the angular momentum is positive, L > 0, then 
because of the modulus signs in Eq. (^) the energy is 
minimised if one takes only positive integers for m. Then 
for < I77I <C 1 degenerate perturbation theory should 
be used with a basis set consisting of all states {Nm}, 

such that L = X]m=0 "^^m- 

The single particle states with n-^ ~ and m > are 



of the form z™exp(— |z| /2) where z 



iy and m 



is the angular momentum quantum number. With the 
use of a computer algebra package we have been able to 
study systems of up to six particles comprehensively. In 
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addition we can prove that the form of the ground state 
holds for an arbitrary number of particles. 

We find that in all cases the groundstate for the at- 
tractive interaction (rj < 0) is = ^^PJ2iLii~\^i\'^) 
where Zc = X]i=i the centre of mass and L is the 
total angular momentum in the system. The contact in- 
teraction energy contribution for the groundstate is in- 
dependent of L, €{N,L) cx r]N{N - l)/2. To prove the 
form of the groundstate wavefunction we show that tpz^ 
is the unique eigenfunction of y = X]i<j ^i^i ~ ^j) 
responding to its largest eigenvalue, Amax- First we note 
that Tpzc is trivially an eigenfunction of V when L — 
since it is the only state in the L = subspace. Since 
the centre of mass coordinate, z^, can be separated out in 
the Hamiltonian, it follows that ipz^ is an eigenfunction 
for any L. 

Let us now work in the basis 



|mi,m2 . . -mj^) 
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where '^rm — L. The matrix elements (m|y|r 
non-negative, and are positive when rrii + rrij = n 
for some i, j, fc, Z, and rup = m'^ for the remaining labels. 
The coefficients {m\'ipz^) are all positive, from which it 
follows that the eigenvector -ipz^ belongs to the largest 
eigenvalue \max- To see this note that Amax can be de- 
rived from variational principle 
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If we take an eigenvector of Amax, and replace all its com- 
ponents by their absolute values, the variational func- 
tional in Eq. (Q) cannot decrease, and so must remain 
at Amax- It follows that there is an eigenvector of Amax 
whose components are non-negative; ipz^ has non-zero 
overlap with this eigenvector, and so must belong to 

Amax ■ 

To prove non-degeneracy of the eigenspace of Amax we 
note that the matrix Vij is "connected" in the following 
way: if we take a basis vector \i) and consider all \j) 
with Vij > 0, then consider all |fc) with Vjk > and so 
on, this includes all basis vectors. If the eigenspace of 
Amax is degenerate then there must be an eigenvector of 
Amax whose components are of both signs (since we can 
choose this eigenvector to be orthogonal to tpz^), and all 
non-zero (if some components are zero, simply add on a 
very small amount of ipz^)- The vector made by taking 
the absolute value of the latter's components will also be 
an eigenvector. The difference in value of the variational 
functional for the two vectors can only be zero if the i- 
th and j-th components have same sign when Vij > 0. 
By connectedness we see that this means all components 
must have the same sign, and hence the eigenspace is 
non-degenerate. This completes the proof that tpza is the 
non-degenerate ground state of the Hamiltonian Eq. (|l|) 
with attractive interaction. 



To determine the degree of condensation, if any, the 
single particle density matrix, p{z, z *) is required for the 
ground state. Yang |2^] showed that off-diagonal long- 
range order is associated with the largest eigenvalue of 
the density matrix (the magnitude of the eigenvalue is 
the fraction condensed) with the 'condensate wavefunc- 
tion' being the associated eigenvector. The notion of off- 
diagonal long-range order is not of such great use for 
trapped atoms, but this definition of the condensate is 
useful in an inhomogeneous setting. The single particle 
density matrix has the form: 



N 

p{z,z'*)^- / I]dz.dz; V( 
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Z,Z2,...ZAr)V'*(z',Z2,..Z7v), 



(5) 



where Q is the normalisation. On integrating we find, 

(iV-l)^-™L! 
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Thus the resulting eigenfunctions and eigenvalues for a 
given m are 



= P. 2 
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If we now consider the case of L = Nq (which in a 
conventional system, e.g. ^He would correspond to q vor- 
tices) then if a condensate exists its eigenvalue will cor- 
respond to m = q. Simplifying we find 



= (1 - 1/N) 



{qN-q)_ 



(Nqy. 
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which can be rewritten as a Poisson distribution in 
limit that iV — > oo. On taking the further limit of g ^ cxd 
the maximal eigenvalue becomes Pq ~ 1/ y/2nq. However, 
the eigenvalues of significant weight are distributed over 
q — yfq "yn < q ~\- ^Jq. This is clearly not the pro- 
nounced peak required for a condensate and is reminis- 
cent of Nozieres and Saint James' fragmented condensate 
pO| . Lest this be thought to be misleading for small 
we note the following results for q = \. We find that 
the eigenvalue where the putative 'condensate' would be, 
pi(q = 1) = e~^, that po(9 = 1) = as well and that 
^2(9 = 1) = 56^^- The 'condensate' is not singled out as 
having a uniquely large eigenvalue. 

The first excited state, is also of interest, as we find 
a rudimentary 'vortex'. The general form is: 

N 

^ = z^-'^^{zi- Zjf with ei (X iV(iV- 2)/2. (9) 

i>j 

Again from symmetry considerations we find that we re- 
quire a minimum of 2 quanta of angular momentum in 
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order to produce an excited state. For L=2, we find that 
there are two possible states for all N and these cor- 
respond to the groundstate and excited state we have 
described above. Because of the separation of the cen- 
tre of mass variables mentioned above, this is in fact a 
general result: for ip^a = then an excited state is $ 
(although we have not proved that it is always the first 
excited state). 

To determine whether the results from the contact in- 
teraction model are likely to be generic or are artifacts 
(for instance of degenerate perturbation theory) , we turn 
to the second model. The Hamiltonian |15|,|l^ (first dis- 
cussed in the context of nuclear physics ^M) describes N 
bosons with attractive harmonic coupling A > 0) (i labels 
the particles): 

^--|^Ev?4^-^ + iB— (10) 



where we enforce the symmetry of the wavefunctions 
at the end of the calculation. In dimensionless units, 
A = A/k, y = {h^/{mk)y/^x and V = {mk/h^y/^V 



N N . N 



which upon rearrangement leads to 



1 ^ 



(1 + A^A) 
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(12) 



Here we note that the problem in d dimensions separates 
into d one-dimensional problems. Hence we will now re- 
strict ourselves to one-dimension for clarity. 

To determine the degree of condensation we again need 
to calculate the single particle density matrix. To do this 
we change variables to the centre of mass coordinate, 

Cn = ^/VNY.i' Vi^ and Ci, {i = 1, N-1) which are ar- 
bitrary but orthogonal to C^r . The Hamiltonian in these 
variables is 



(1 



^Ecf 



NA 



Cat I 



(13) 



This leads to the groundstate wavefunction having the 
form 



-i(l+JVA)i/^C%-^(l-(l+A'A)'/')C^ 



(14) 



where C = {Cii --Xn}- The corresponding frequencies 
are ei = {1 + NAy/'^/2 ior i ^ N and cn = 1/2. We 
now consider the model in d > 2 with attractive interac- 
tions in the fixed angular momentum subspaces. In the 
groundstate we find that all angular momentum is ab- 
sorbed by the centre of mass variable. This may be seen 



by noting that the centre of mass oscillators (associated 
with the different components of the centre of mass mo- 
tion) have a lower associated frequency than the other 
oscillators describing relative motion. The physical in- 
terpretation is straightforward: relative motion requires 
more energy as work must be performed against the at- 
tractive interactions. (The converse will hold true for 
repulsive interactions). Hence, for two-dimensions and L 
quanta of angular momentum we can immediately write 
down the wavefunction, 



^ = ^Lg-i(l + iVA)^''"N|%-i(l-(l + iVA)^/")|^c 



(15) 



In the thermodynamic limit it can be shown that the 
contribution to the single particle density matrix of the 
exponential term associated with Zc is negligible. Hence, 
surprisingly, the density matrix reduces to that of the 
contact interaction model. The groundstate wavefunc- 
tions will therefore be the same, as will the properties of 
the single particle density matrix. 

We shall now show that these systems are condensed 
(at least under some conditions) when the interactions 
are repulsive. Thus the lack of condensation is not due 
to peculiarities of the models in general, but of the at- 
tractive interactions in particular. 

The repulsive case has also been considered, and will be 
discussed elsewhere. We find no vortex lattice, although 
a Laughlin like state exists for some special values of L. 

Firstly consider the contact interaction model when 
there are N quanta of angular momentum in the system. 
Conventionally there would then be one vortex in the sys- 
tem, and one might expect that the ground state (subject 
to the constraint of the amount of angular momentum) 
would be: 



N 



N 

n 
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(16) 



We conjecture the following form for the CIM: 

N 

rL=N = n ([^1 + • • • + ^« - Nz,]e-\^-\'/') (17) 

We have demonstrated that this form is correct by ex- 
plicit calculation on systems of up to 6 bosons. The 
physical interpretation of this wavefunction is that the 
bosons are rotating around the centre of mass, which 
would be a condensate if the centre of mass were a c- 
number. We will now show that in the thermodynamic 
limit the corrections to full condensation are 0{1/N). 
Consider the density matrix constructed from the wave- 
function Eq. (irr 
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X - ^ -z \ \ uj — -z 



N 



n 



X I M |c^ + ^-z,j + le 
i=2 



where uj = l/-^X]iI=2 -Zj-Then to separate the integration 
over the different Zi we introduce a delta function for the 
centre of mass variable: 



(18) 



1—1 duJxdiOy 5 ojx 
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and use the integral representation 
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Upon substitution and integration we find that, in the 
limit that N oo, to accuracy 0{1/N'^): 



p(2,z'*)=^^(z)p„^*(z') 



2 ■ 
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(21) 



Thus there is a condensate with eigenvalue 1 — 2/N, in the 
state 2, which is fully condensed in the thermodynamic 
limit. The corrections of 1/A^ are in the states 1 and z^. 
In addition a Laughlin state is empirically found in small 
systems to be the ground state for L = N{N — 1). 

Turning to the harmonic interactions model, we note 
that the centre of mass oscillator has a higher frequency 
than the others when the interaction is repulsive. In that 
case the other oscillators will be populated in preference 
when minimising the energy subject to L (a multiple of 
N) quanta of angular momentum. Now, the other oscil- 
lators are degenerate and the centre of mass factor is not 
being multiplied by a large multiple of the centre of mass 
coordinate, which it was in the attractive case. The lat- 
ter implies that the centre of mass factor is irrelevant in 
the thermodynamic limit. Thus the ground state reduces 
to a single particle form [Eq. (||)] and hence the answer 
will be a condensate into the state z^^^. 

In conclusion we have shown that for a model with 
degenerate ground states in the absence of interaction, 
there is no condensate formed when weak interactions 
are incorporated. Consequently, in this particular case, 
there is no vortex lattice. In a different model, which 



does not have degenerate ground states, we have shown 
that the particles are uncondensed when given an exten- 
sive quantity of angular momentum. In both cases the 
angular momentum of the system resides in the centre 
of mass motion, in contrast to the more familiar case of 
repulsive interactions. This leads to the general hypoth- 
esis: attractive bosons do not condense in the presence of 
single particle degeneracy and their angular momentum 
resides in the centre of mass motion. The investigation of 
rotating "^hi might be fruitful in exposing an uncondensed 
'ground state'. Repulsive interactions in the same models 
lead to condensed ground states, so showing that attrac- 
tive interactions do indeed lead to different physics. 
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